Kink dynamics, sinh-Gordon solitons and strings in AdS3 
from the Gross-Neveu model 
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Guided by a study of kink-antikink scattering in the Gross-Neveu model and other known so- 
lutions of the Hartree-Fock approach of a particularly simple type, we demonstrate a quantitative 
relationship between three different problems: Quantized 1+1-dimensional fermions in the large N 
limit, solitons of the classical sinh-Gordon equation and classical strings moving in 3-dimensional 
anti de Sitter space. Aside from throwing light on the relationship between quantum field theory 
and classical physics, this points to the full solvability of the dynamical iV-kink-antikink problem in 
the Gross-Neveu model. 

PACS numbers: 11.10.Kk,11.10.St,11.25.Pm 



I. INTRODUCTION 

The simplest variant of the Gross-Neveu (GN) model 
family Q consists of N species of massless Dirac fermions 
in 1+1 dimensions, interacting via a scalar-scalar four- 
fermion interaction, 



iV 



k=l 



N 



(1) 



This model has asymptotic freedom, no scale, and a dis- 
crete chiral symmetry tp — > 75^ which gets broken spon- 
taneously in the vacuum, yielding a dynamical fermion 
mass m via dimensional transmutation. Throughout this 
paper we will be exclusively dealing with the 't Hooft 



limit TV — > 00, Ng = const. 



As is well known, the 



attractive interaction gives rise to a marginally bound 
scalar fermion- ant ifermion state (the a-meson) with mass 
2m and to a rich variety of multi-fermion bound states 
(baryons) [!, H|. Moreover, the model features a non- 
trivial phase diagram as a function of temperature and 
chemical potential with three distinct phases (massless 
and massive Fermi gas, baryon crystal) meeting at a tri- 
critical point [f| . It is of interest not only as a toy model 
for strong interaction particle physics, but also because 
of its almost literal recurrence in condensed matter sys- 
tems such as conducting polymers, carbon nanotubes or 
quasi-one-dimensional superconductors @. In the large 
N limit, baryons as well as baryonic matter and the phase 
diagram can be determined with scmiclassical methods, 
notably the relativistic Hartree-Fock (HF) method 0. 
Relativity enters in two ways — use of the Dirac equa- 
tion instead of the Schrodinger equation, and taking into 
account the filled, interacting Dirac sea. Recently we 
have started to address time-dependent questions by gen- 
eralizing this approach to time-dependent Hartree-Fock 
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(TDHF) In a first step, the boosted baryon was con- 
sidered, demonstrating a covariant energy-momentum re- 
lation and deriving exact structure functions as fermion 
momentum distributions in the infinite momentum frame 
. In the present work, we proceed to the next level of 
complication and study baryon-baryon scattering, an is- 
sue which has not yet been discussed in any detail in 
the large N limit. This may well be the first concrete 
realization of Witten's vision about baryon-baryon inter- 
actions in the large N limit, originally developed in a 
non- relativistic context [lOj . Aside from the genuine in- 
terest in the scattering problem of composite, relativistic 
bound states, we hope to get further insights into the 
mathematical structure of the GN model by enlarging 
the spectrum of questions addressed in the framework of 
the TDHF approach. 

Before embarking on this problem, it may be worth- 
while to recapitulate some of the experience gained with 
previous applications of the HF approach. The basic 
mathematical problem can be formulated in a single line 



(i£ - S) V Q = 0, 



s 







(2) 



where the ip a are c-number spinors (single particle wave 
functions). Since the sum over occupied states includes 
the Dirac sea, this is an infinite system of coupled, non- 
linear partial differential equations. At finite tempera- 
ture and chemical potential, a similar formula applies, 
but the sum includes thermal occupation numbers. Ex- 
cept in the case of a homogeneous mean field S which 
acts like a dynamical mass, the solution of Eqs. ((2J is 
highly non-trivial. Nevertheless, for the model defined in 
Eq. ([1]), closed analytical solutions have been found in 
all cases studied so far. If one examines in detail how 
self-consistency is achieved, one notices that all known 
solutions fall into two classes: Type I solutions are such 
that the contribution from every occupied single particle 
state is either zero or proportional to the full mean field 
S, 



4>a4>a = KS. 



(3) 
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The self-consistency condition then becomes space-time 
independent, 



Ng 2 J2^a = 1. 



(4) 



Aside from the trivial case of the vacuum, type I solutions 
have been found for the kink at rest 0, [ill or in flight Q 
and the kink crystal at zero temperature [l2[ relevant for 
the ground state at finite density. In type II solutions, 
the scalar condensate for each occupied state involves two 
different space-time functions. Without loss of generality, 
this can be written as 



V> Q -0a = KS + X' a S', 



(5) 



where S is the self-consistent potential and S' a second, 
different function. This yields two space-time indepen- 
dent self-consistency conditions, 



(6) 



Known type II solutions are the kink-antikink baryons 
[E EH, the kink crystal at finite temperature [l3| and 
the time-dependent breather Q to be discussed below in 
more detail. Incidentally, all HF solutions of the massive 
GN model (i.e., Lagrangian (fTJ) supplemented by a term 

—mo Y^k=l ^ki^k) are also of type II [UHIU, and no solu- 
tion of type III or higher is known in the GN model which 
would require more than two functions and hence more 
than two space-time independent self-consistency condi- 
tions. As will become clear later on, this classification of 
HF solutions is useful if one wants to relate the quantum 
theory to classical GN models with a small number of 
flavors. 

In this paper, we shall focus on the dynamics of kinks, 
supplementing the known single kink and kink crystal so- 
lutions by kink-antikink scattering and briefly comment- 
ing on the generalization to N kinks and antikinks. It 
turns out that this restriction to kink dynamics is at the 
same time a restriction to type I HF solutions. For this 
class of particularly simple (though exact) solutions of 
the large N, massless GN model, we will identify an ef- 
fective bosonic theory whose soliton solutions are closely 
related to self-consistent HF potentials, using methods 
developed by Neveu and Papanicolaou in their proof of 
integrability of the classical one- and two-flavor GN mod- 
els [17| . The relevant equation is the sinh-Gordon equa- 
tion, the hyperbolic version of the more familiar sine- 
Gordon equation. This in turn is the key for mapping 
type I solutions of the GN model onto classical string 
theory in 3-dimensional anti de Sitter space (AdS.-s), fol- 
lowing recent work of Jevicki and collaborators [18|-|20|. 
In this way we hope to show that in spite of the simplic- 
ity of the GN model, its mature age and the considerable 
amount of work devoted to it, there is still room for new 
insights and surprises. 



This paper is organized as follows. In Sec. UH we solve 
the kink-antikink scattering problem with the help of the 
known breather. Sec. IIIII aims at identifying the sinh- 
Gordon equation as underlying effective bosonic theory. 
In Sec. IIV1 we show how to map solutions of the GN 
model onto strings moving in AdS,3 and illustrate this 
mapping with a few examples. In Sec. [V] we summarize 
our findings and identify promising directions for future 
work. 



II. KINK-ANTIKINK SCATTERING 

Our starting point is the kink-antikink breather so- 
lution of the GN model discovered by Dashen, Hass- 
lacher and Neveu (DHN) |[. Since the inverse scattering 
method is not available for time-dependent semi-classical 
solutions, these authors had to guess the scalar potential 
by analogy with the well-known sine-Gordon breather. 
They then show that a self-consistent solution can be 
found whose quantized fluctuations yield the spectrum 
of kink-antikink type baryons. At the end of their paper, 
DHN mention that the imaginary choice e = i/v of a cer- 
tain parameter governing the frequency of the breather 
should describe scattering of a kink-antikink pair with ve- 
locities ±w in their center-of-mass (cm) frame. We take 
up this suggestion here and analyze the scattering of a 
kink and an antikink in greater depth, using the language 
of the TDHF approach. Before doing this however, let 
us review the results for the DHN breather which are 
relevant for the analytic continuation in e. This is also 
necessary because the formulae given in 0] do not seem 
to be fully consistent, either due to typos or to some 
unspecified conventions. 

The breather is a solution of the Dirac equation 



(i@ -S)ip = 
and will be given in the representation 



(7) 



7 =-CTi, 



7 1 = if7 3, 75 = 7°7 1 = -°2 (8) 



of the 7 matrices. Units in which the vacuum fermion 
mass is 1 will be used throughout this work. The scalar 
potential can be written as 



with 



S = l+£f 2 + Vf4 



h = f4CosQt, 

,fi = (cosh Kx + a cos fit + b) . 



(9) 



(10) 



The parameters appearing here arc defined and related 
as follows, 



£ = -2a, T) 



K = eCl, 



-bK 2 



(11) 
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a is the solution of the equation 



= rf 



1-b 2 



n 2 



(12) 



where we choose the positive square root. S is fully speci- 
fied if the parameters e and b are given. The Dirac equa- 
tion ([7]) with potential © has two types of solutions. 
First, there is a continuum of wave- like solutions of the 
form 



ipk 



<7v 



i(kx— ut) 



(13) 



with ui = zt^/k 2 + 1. Like in the single baryon problem, 
these scattering states have no reflected wave, showing 
that also the time-dependent potential is transparent. 
DHN give the explicit form of ipk in terms of /2, f± from 
Eq. (TIU)) and two further functions 



fx = f 4 sinh Kx, f 3 = f 4 sin M. 



(14) 



We have redetermined the coefficients with the following 
results, 



Uk 



<7v 



A4 
A4 



iK 



2k 
1 + ik 



ia ^ ail 1 + 
k 2k uj 



ik it] 
~ h+ 2k U 



K i 



2k 



a i — i 
k ui 



■h 



aO , r) i - k 
"2fc /3 2k uj U 



(15) 



They deviate from Eq. (4.6) of Ref. ||, but one can easily 
check that Eqs. (|13|15|) do solve the Dirac equation. The 
second type of solution are bound states. DHN find two 
distinct states which can be expressed in terms of the 
functions 

<Pi = fi cos &t /2 cosh K x/ 2, 
¥2 = fi sin Ot/2 sinh Kx/2, 
<P3 = f4smttt/2 cosh Kx/2, 
(fi = fi cos Qt/2 sinh Kx/2. (16) 

Two orthogonal solutions of the Dirac equation are 



(i) 



4 2) 



A/o 



C2f2 



C4V4 
C 3 (f 3 



(17) 



Here, our coefficients 



C2 = 



C3 



C4 



l-b-n/2 
a — rj/2 
iQ 1 + b — a 
Tl + b + rj/2' 
Kl-a-b 



77/2 



(18) 



differ from those of DHN in the overall signs of c 3 and 
C4. Self-consistency of this solution can be established as 
follows: DHN find that ipktpk for each occupied contin- 
uum state yields two contributions, one proportional to 

5 and one proportional to / 4 . The discrete states yield 
only a single contribution proportional to f 4 . By relating 
the parameter b to the occupation fraction of the discrete 
states, one gets a self-consistent result where the f 4 terms 
cancel. According to the classification of TDHF solutions 
given in the introduction, this shows that the breather is 
a type II solution of the GN model. Only for the value 

6 = 0, the contributions proportional to f 4 would vanish 
and the breather would be a type I solution. However 
b = is ruled out since it is incompatible with a real 
potential S. 

We now turn to the kink-antikink scattering problem. 
Following the suggestion of DHN, we analytically con- 
tinue the above breather solution to the value e = i/v. 
This gives rise to the following changes in our formulae 
for the potential and the spinors, 



K 

n 

cos ttt 
sin fit 



VT^u 2 
2iv 



VT^j 2 ' 

2vt 



cosh ■ 



-i sinh ■ 



2vt 



— * 

1* 



b 2 



l-v 2 



(19) 



The continuum states now have to be normalized to free 
spinors at t — > — 00 using 



2k 2 



4fc 2 + K 2 ' 



(20) 



The discrete states are square integrable and normalized 
to 1 by 



A/o 



K(l + a + b) 



(21) 



Let us now consider the issue of self-consistency for the 
scattering problem. The contribution from the negative 
energy (uj < 0) continuum to the scalar potential is found 
to be 



= SNg 2 



A / 2 dk 



1 



A/2 27T y/k 2 



with 



h(k) 



+f4Ng 2 



2bv 2 



dk 
2^ 



v / r+fc2(i - v 2 )(i + k 2 - v 2 k 2 ) 



h{k) (22) 



(23) 



4 



The first term by itself gives the self-consistent result 
owing to the vacuum gap equation, 



1 = Ng 2 



A / 2 dk 



1 



-A/2 27T y/k~2~+T' 



(24) 



This observation is the same as for the breather. Differ- 
ences arise in the treatment of the discrete states. Since 
there are two orthogonal bound states, there is an am- 
biguity how to fill these states in the TDHF approach. 
DHN argue that in the case of the breather one should 
take the linear combinations t/j^^ip^ and fill them with 
N and uq fcrmions, respectively. We have to reconsider 
this prescription in the light of the scattering problem. 
Appropriate initial conditions for the TDHF equation are 
an incoming kink and antikink with prescribed baryon 
numbers. In order to prepare such an initial state, we 
need to know which linear combinations of the two bound 
states goes over into the kink or antikink (zero-energy) 
valence state for t — > — oo. Likewise, in order to be able 
to interpret the final state of the scattering process, we 
need the corresponding analysis at t — > oo. Consider an 
arbitrary (normalized) linear combination of tp^i^^t 



(1) 



m4 2) 



(25) 



In order to exhibit the asymptotic behavior of this wave 
function for x — > ±oo,i — > ±oo, we introduce labels for 
the incoming/outgoing wave from/to the left/right as fol- 
lows, 



il 

ir 
ol 
or 



— OO, X 
— OO, X 

+oo, X 
+oo, X 



-00 

-oo 
-oo 

-00 



(26) 



Using this shorthand notation, we find the asymptotic 
expressions (omitting a common normalization factor 



■0O.il 
00,11 
0O,ol 
V>0,or 



1 



A(l - c 4 ) + ifi(c 2 + c 3 ) 
coshz + V -«A(c 2 - c 3 ) + /i(l + c 4 ) 

1 ( A(l + c 4 ) - i//(c 2 - c 3 ) 
coshy_ \ «A(c 2 + c 3 ) + ^(1 - c 4 ) 

( A(l - c 4 ) - i[i(c 2 + c 3 ) 



coshy + \ «A(c 2 - c 3 ) + + c 4 ) 

( A(l + c 4 ) + in(c 2 - c 3 ) , 
-iA(c 2 + c 3 ) + /i(l-c 4 ) ' AZ<) 



with 



cosh ; 



V± 



- (K(x + vt) ±lna) . 



z± = — (K(x — vt) ± In a) 



(28) 



We introduce two orthonormal discrete states with pa- 
rameters (Ai, jxi) and (A 2 , ^2) and occupation N\, In 



order to match the initial conditions of the kink-antikink 
scattering problem, we require that state 1 has only in- 
cident fcrmions from the left and state 2 only incident 
fcrmions from the right. Eqs. (|27[) then imply the two 
pairs of homogeneous, linear equations 

= Ai(l + c 4 ) -i/x(c 2 -c 3 ) 

= iAi(c 2 + c 3 ) -c 4 ) 

= A 2 (l - c 4 ) +i/x 2 (c 2 + c 3 ) 

= -iA 2 (c 2 -c 3 )+M2(l + c 4 ) (29) 

The conditions for the existence of a non-trivial solution, 



det 



det 



1 + c 4 -i(c 2 - c 3 ) 
«(c 2 + c 3 ) 1 - c 4 

1 - c 4 i(c 2 + c 3 ) 
-i(c 2 - c 3 ) 1 + c 4 



(30) 



are indeed satisfied for arbitrary v, b. We can then deter- 
mine the ratios /i;/Ai, 



Mi 
Ai 

m 

A 2 



. I + C4 
-i . 

c 2 - c 3 

1 — c 4 

I . 

c 2 + c 3 



(31) 



The right-hand side is purely imaginary. This implies a 
vanishing contribution to the chiral condensate from the 
discrete states which is proportional to A*/i + /i*A, as 
can easily be checked. Self-consistency in the scattering 
problem is therefore only possible for b = 0. In this 
particular case, the coefficients Cj are 



c 2 = 1, c 3 = -c 4 



1-a 



so that Eqs. (|3 1 [) simplify to 

Mi _ _• _ _M2 
Ai A 2 



(32) 



(33) 



A consistent choice of parameters is therefore b = and 



Ai — A 2 



1. 



/'1 



"M2 



(34) 



For these parameters, one finds that state 1 has only an 
incoming wave from the left and an outgoing wave to 
the right, state 2 only an incoming wave from the right 
and an outgoing wave to the left. The valence fermions 
are exchanged between the two scatterers, presumably a 
consequence of the fact that the potential is transpar- 
ent. These valence fcrmions do not play any role in the 
self-consistency issue. We can prescribe the number of 
valence quarks independently for kink and antikink as 
Ni,N 2 where < Ni < N. In the final state, iVi, iV 2 
are then simply exchanged. We recall that the kink (or 
antikink) with Ni valence quarks carries fermion number 
Nj = Ni — N/2 as a result of induced fermion num- 
ber in a topologically nontrivial background potential 
[2lL [22l | . Expressed in terms of reduced fermion number 



5 



Nj/N = Ni/N - 1/2 = ui 
thus reads 



1/2, the scattering process 



K(v 2 -l/2) + K(v 1 
The time delay 



1/2) -►A'Ci/i -1/2) 



A* 



lnu / — 

— VI 

v 



1/2). 

(35) 



(36) 



is independent of the fcrmion numbers and negative, indi- 
cating a repulsive interaction between kink and antikink. 
Note that the value 6 = implies that the kink-antikink 
scattering solution is of type I and hence significantly 
simpler than the breather. This is also reflected in the 
form of the scalar potential, 



S 



v cosh Kx — cosh Kvt 
v cosh Kx + cosh Kvt ' 



K 



(37) 



Incidentally, Eq. ([37]) agrees with a result quoted in 
Ref. [13] without derivation. The fermion density for the 
discrete states is given by 



Pl,2 



vK [v + cosh K(x ± vt)] 
2(v cosh Kx + cosh Kvt) 2 



(38) 



where the + sign belongs to state 2 moving to the left, 
the — sign to state 1 moving to the right, cf. Eq. ([Ml) . 
Expressions ([3"5|) have been normalized to 1 and have to 
be multiplied by the occupation numbers to get the true 
fermion density. The density ipt^k of the negative energy 
continuum states is 



Pk 



1 



(k 2 v 2 



2v (v + cosh Kx cosh Kvt) 
- k 2 — 1 ) (v cosh Kx + cosh Kvt)' 
2kv 2 sinh Kx sinh Kvt 



(39) 



\fl + k 2 (k 2 v 2 — k 2 — 1) (v coshifa: + coshKvt) 

Integrating over k and subtracting the vacuum contribu- 
tion yields the simple, finite result 

dk 1 . 

7r(Pk ~ !) = -7;(Pi + P2)- 

The total fermion density per flavor is thus given by 



PS=\v\ 



Pi 



1 



V2 - 



P2- 



(40) 



Integration over x finally yields the (conserved) total 
fermion number, 



N f f°° 

— = I dxp f (x 7 1) = Vi + v 2 - 1 

^ J — oo 



(41) 



The kink-antikink scattering process is illustrated in 
Figs. HHSl Fig. [T] shows, by means of the scalar potential 
S, that kink and antikink approach each other, are re- 
pelled and bounce back. This picture is independent of 
the baryon numbers involved. Figs. [2] and [3] exhibit the 
baryon density for baryon-antibaryon and baryon-baryon 
scattering. In the first case one clearly sees the exchange 
of fermions during the collision. 




FIG. 1: Scalar potential for kink-antikink scattering at v 
0.5, showing a repulsive interaction. 




FIG. 2: Fermion density for kink-antikink scattering. Param- 
eters: v — 0.5, V\ = 1, v% = (baryon-antibaryon collision). 



III. EFFECTIVE BOSONIC FIELD THEORY 

A big hurdle in applying the TDHF approach is to find 
the self-consistent potential. Obviously it would be ex- 
tremely nice if one could write down a closed equation 
satisfied by the potential. This has been achieved ap- 
proximately in the case of the massive chiral GN model 
near the chiral limit, using the derivative expansion [l4j |. 
Such classical field equations may be regarded as equa- 
tions of motion of an effective scalar theory where the 
fermions have been "integrated out" , in the language of 
the path integral. The derivative expansion is of little 
use in the discrete chiral GN model where the baryons 
remain localized in the chiral limit. Here we address the 
question whether one can nevertheless identify an effec- 
tive bosonic theory for the non-chiral GN model. 
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FIG. 3: Fermion density for kink-antikink scattering. Param- 
eters: v — 0.5, V\ = 1, i/2 = 1 (baryon-baryon collision). 



Since the tanh kink appears in 4 theory, one may be 
tempted to postulate the equation 

d^S + 2S 3 - 25 = 0. (42) 

It is indeed solved by the boosted kink tanh(7(a; — vt)). 
However, this cannot be correct because Eq. (|42[) has no 
solitons in the strict sense of the word: If one scatters 
a kink and an antikink, a numerical study shows that 
they interact inclastically and do not keep their tanh 
shape [HI • A simple bosonic theory which does not have 
this problem can be found for type I solutions only, to 
which we restrict ourselves in this section. These solu- 
tions are specific for kink dynamics (kink, kink crystal, 
kink-antikink scattering), and therefore only exist in the 
masslcss GN model. For type I solutions we can write 

S = l a ^a (43) 

for every occupied state a with nonvanishing ipatpa- If we 
pick any such state, the TDHF problem evidently reduces 
to the N = 1 classical GN model, i.e., a non- linear Dirac 
equation 

(0 - tatpaTpa) 0a = (44) 

with c-numbcr spinors. We suppress the label a from 
now on. Ncvcu and Papanicolaou [l?} have proven the 
integrability of the classical N = 1, 2 GN models long 
ago. The part of their work dealing with N = 1 in fact 
contains the answer to the question about the effective 
bosonic theory. To demonstrate this fact we closely follow 
their work in the present section. 

To this end it is useful to switch to a different repre- 
sentation of the 7-matrices in which the upper and lower 
components of the Dirac spinor ipi , ip2 have definite chi- 
rality, 

7° = o-ij 7 1 = i( ?2, 75 = -0-3 ■ (45) 



Using light-cone coordinates 

z = x — t, z = x + t (46) 

and the comma-notation for partial derivatives, the non- 
linear Dirac equation assumes the simple form 

-2iip 1>z = S^ 2 , 

2^2,2 = S^i, (47) 

where 

S = l{^2+^i). (48) 

The Dirac equation expresses ipi^ z and xp2,z m terms of 
V'l) V^- What about the other derivatives, ijii s and ^2,z? 
Following Neveu and Papanicolaou, one first derives the 
identities 

Styi,z - S^ipi = -ihi£ip2, 
Sip 2 , z -S tZ ip 2 = -ih 2 tyi, (49) 

with 

hi = i ^1,2 - V'iV'm) i 

h 2 = i tyW2,z ~ M*2, z ) ■ (50) 

Since the right-hand side depends again on the unknown 
derivatives ipi.zT^.z and their complex conjugates, it 
looks as if the goal of expressing these derivatives through 
■01, -02 had not been reached. However, simple algebra 
shows that 

h hz = = h2,g (51) 

so that hi t 2 can only depend on either z or z. If these 
functions would not be constant, they would describe 
fermion bilinears propagating at the speed of light. In 
the absence of massless particles in the model we con- 
clude that hi, are constant for physically sensible so- 
lutions. This enables us to express all derivatives of ip in 
terms of ip itself as follows, 

ip,g = Cii>, i>, z = C 2 0, (52) 

with the matrices 



\ -iS/2 ) ' 




The integrability condition for the system (|5"2"|) reads 

Ci, z -C 2 ,g + [Ci,C2]=0 (54) 

or equivalently 

SS, zz -S, z S tZ -^S 4 = hih 2 £ 2 (55) 
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As the right-hand side is constant, we have succeeded 
in deriving a non-linear partial differential equation for 
the self-consistent TDHF potential, at least for type I 
solutions. If we disregard the kink crystal for the mo- 
ment and focus on localized solutions where S reaches 
its vacuum value asymptotically, we can even determine 
the constant. In units where S = m = 1 in the vacuum, 
we infer from Eq. (|55[) that 



hih 2 £ 2 = -- 



(56) 



In normal coordinates, the final equation for S reads 

Sd^S ~ d^Sd^S + S 4 - 1 = 0. (57) 

It is satisfied by the kink and the kink-antikink solutions 
in arbitrary Lorentz frames, but also by the kink crystal 
for which we have not derived the value of the constant, 
Eq. (|56[) . One can write down a simple (although singu- 
lar) Lagrangian which yields Eq. ((57)) as Eulcr-Lagrangc 
equation, namely 



1 



C =g2 {d^Sd^S -S 4 -l) 



(58) 



However, this should not be interpreted as effective 
bosonic field theory for the GN model. If one derives 
the Hamiltonian density H from C in the usual way, one 
does not get the correct energy density. The reason is 
presumably the fact that we already made use of proper- 
ties of the solution when deriving Eq. (f5"Tj) . notably the 
fact that we are dealing with a type I solution of the 
TDHF problem. In this way we are not really able to in- 
tegrate out the fermions in full generality and construct 
the effective bosonic action for arbitrary scalar fields. 

Eq. (|57|) does not yet resemble any of the well-known 
equations with solitonic solutions. The closest we could 
come to a more familiar looking form was by means of 
the change of variables 



17] 



5- 



= lnS 2 



(59) 



which reduces Eq. (|57|) to the sinh-Gordon equation, 

<9 M <9' 1 + 4sinh0 = O. (60) 

However we loose the information about the sign of S 
and hence also the Z2 chiral symmetry in this nonlinear 
transformation - the two vacua S = ±1 are mapped onto 
the same value 9 = 0. Since the zero's of S give rise to 
singularities of 0, it is easy to reconstruct a solution S of 
Eq. (|57|) from a singular solitonic solution 8 of the sinh- 
Gordon equation, so that the mapping is nevertheless 
quite useful. With this caveat, the matrices C±,C2 and 
the linear equations (|5"2"j) can then be identified with the 
Lax pair of the sinh-Gordon equation. 

Notice that the coefficient 4 in Eq. (|60p has a sim- 
ple physical interpretation: The linearized sinh-Gordon 
equation, 



yields the Klein- Gordon equation for a scalar meson with 
mass 2 in units of the fcrmion mass. This can be identi- 
fied with the a meson of the massless GN model. The re- 
lation between kink, a meson and sinh-Gordon equation 
in the massless GN model is analogous to the relation be- 
tween the light baryon, the 7r meson and the sine-Gordon 
equation in the massive 2d Nambu-Jona-Lasinio model 
(NJL2) close to chiral limit (derivative expansion 
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In the latter case, this was interpreted as the simplest re- 
alization of the Skyrme model in the sense that baryons 
are topologically non-trivial excitations of the pion field 
[24j| . In the GN model, we have now identified a similar 
picture for the case of a discrete chiral symmetry, the 
baryon emerging as a large amplitude excitation of the a 
field. The similarity between the baryons in both cases 
is particularly striking if one writes the sine-Gordon and 
sinh-Gordon solitons in the following form: 

1. Massless GN model (exact, m a = 2), 



S 2 = e\ 

= d^e 



= —4 artanh 
n 2 sinh 9. 



(62) 



2. Massive NJL2 model (leading order derivative ex- 
pansion, TOtt = 2^/7 with 7 the confinement param- 



eter 0), 

S-iP 




4 arctan e r 



m 2 sin ( 



(63) 



We can also write down a common formula for topological 
baryon number, 



In $(00) - In $(—00) = 2iriNB, 



(64) 



with $ = S for the GN model and $ = S - iP for 
the NJL2 model. However, in the GN case, this only 
determines the non-integer part of the induced baryon 
number [25[ (kink and antikink give =f1/2 although they 
have the same value of induced baryon number of —1/2) 
so that the analogy should perhaps not be overrated. 

The kink belongs to the class of "traveling wave solu- 
tions" of the sinh-Gordon equation. Kink-antikink scat- 
tering is an example of a "functional separable solution" 

M 



where 



8(x, t) = 4 artanh [f(t)g(x)] 

(f, t ) 2 = Af + Bf + C, 
-(9,x) 2 = Cg* + (B + 4)g 2 + A. 



(65) 



(66) 



{d^ + 4) e = 0, 



(61) 



The general TV soliton solution is also known for the sinh- 
Gordon equation 0, H3| and is a likely candidate for the 
TDHF solution of the GN model with N kinks and an- 
tikinks. Since the Lax pair including the spinor wave 
functions are known, all what one would have to do is 
find the bound state solutions and verify self-consistency. 
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Prom the point of view of particle physics, this is per- 
haps somewhat academic because it describes a scatter- 
ing problem with N incident baryons. Nevertheless, it 
would be challenging to solve the relativistic 7V-baryon 
problem exactly and analytically in terms of the elemen- 
tary fermion constituents, including the Dirac sea effects. 



IV. RELATION TO STRINGS IN ADS 3 

Jevicki et al. have explored the close relationship be- 
tween sinh- Gordon theory on the one hand and classical 
strings in AdS 3 on the other hand [H-HSl. Since the GN 
model can also be related to the sinh-Gordon model, this 
should enable us to map type I solutions of the GN model 
onto solutions of classical string theory. From the physics 
point of view, a stringy interpretation is not obvious for 
a non-gauge theory like the GN model. It may be of 
interest to see explicitly how such a mapping works. 

Let us briefly review the required string theory back- 
ground following Ref. [lH. The AdS3 target space 
is parametrized by the embedding coordinates Y a (a — 
-1,0,1,2) in R 2 ' 2 , 



Y-Y 



-Y 1 



Y 2 



Y 1 



-1. 



(67) 



The string equation of motion in conformal gauge reads 

(d = d z ,d = dz) 



ddY - (dY ■ dY)Y = 



(68) 



and has to be supplemented by the Virasoro constraints 

(dY) 2 = = (BY) 2 . (69) 

The crucial instrument in relating strings in AdS3 to 
the sinh-Gordon equation is the Pohlmeyer reduction 
[28j . using the factorization of the AdS3 isometry group 
SO(2,2) into SO(2,l)xSO(2,l). The string coordinates 
are expressed through two auxiliary spinors <p, x as fol- 
lows, 



Zi = Y_ 1 +iY Q = 4>lxi - 02X2, 
Z 2 = Y 1 +iY 2 = (f>2X*i - <PiX 2 - 

Here, </> and x are normalized according to 



1 



>l9l 



'J 2 V2 



XiXi - X2X2 



and satisfy the linear system of equations 



X. 



M4, 

BiX, 



4>,z 

X..2 



= A 2 <f>, 
= B 2X , 



(70) 



(71) 



(72) 



where now z, z are light-cone coordinates derived from 
worldsheet space-time coordinates. The matrices depend 
on two arbitrary functions u(z),v(z) which we choose 
as u ~ 2, v = — 2 for the sake of simplicity. For these 
particular values the normalization of the light-cone co- 
ordinates agrees with Eq. (|46[) if we identify x, t with 



the worldsheet coordinates. 
Eq. ([72]) are then given by 



The matrices introduced in 



A 2 
Bo 



—i\c+ ia^ — Ac. 
— ia^ — Ac_ iAc+ 

ic+/\ —ia tZ — C-/\ 
ia^z — c_/A — ic + /\ 

— iAc_ ia^z — Ac + 
— ia^z — Ac + iAc_ 

ic_/A — ia. 2 — c+/A 
ia, z — c+/A — ic-/X 



with 



c± 



-a/2 ± e a/2^ 



(73) 



(74) 



a(z, z) is a function which is expressible in terms of 0, but 
this relation will not be needed here. The integrability 
condition for both systems of linear equations (|T2")) yields 
the sinh-Gordon equation, 



a. 



sinh a. 



(75) 



It is straightforward to verify that a pair of spinors </>, x 
satisfying the normalization conditions (|7ip and the lin- 
ear sytem (|72|) yields a solution of the string equations 
(|67H69p . the link between the two problems being given 
by Eqs. J7DD. 

The way to relate the GN model to strings in AdS3 
is to find a gauge transformation between two different 
representations of the sinh-Gordon Lax pair, Eqs. (|52I53[) 
and Eqs. (|72l73p . where we have to expect complications 
due to the fact that the mapping from the GN model to 
the sinh Gordon model was not one-to-one. Using the 
paramctrization 



S = ±e e ' 2 



(76) 



depending on the sign of S, let us first consider the case 
S > 0. The gauge transformation from (Ci, C 2 , ?/>) to 
(Ai,A 2 ,cj)) is given by 



A l = 
A 2 = 



Q(d -d)Q-\ 

n(c 2 -d)n-\ 



(77) 



A\. 2 depend on a and a scale parameter A, Ci,2 depend 
on S = e B l 2 and the constants hi£, h 2 l related through 
Eq. (|56[) . For future convenience, we introduce a param- 
eter C, (to be identified with the spectral parameter of the 
GN model TDHF solutions) via 



M = 2C 2 , h 2 £ = ~—r 



1 

8C 2 ' 



(78) 



Q, can be found in 2 successive steps: First, make a local 
Abelian gauge transformation to render C\ , C 2 traceless 
(Ai,A 2 are in the SO(2,l) Lie algebra). This is achieved 
by the choice 



ft! =S 



-1/2 



(79) 
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Secondly, perform a global non-Abelian gauge transfor- 
mation 



(80) 



The product 



A- 1 1 
iX- 1 -i 



S~ 1/2 (S > 0) (81) 



then transforms C\ , C2 onto A\ , A2 provided we relate 
the field variables (0, a) and spectral parameters ((, A) 
as follows, 



A 



4C 2 



(82) 



As will be seen later on in concrete examples, the correct 
sign for negative energy continuum states is A = — 2Q. 
Hence the mapping from the GN solution ip to the first 
string a model spinor <f> is 



4> = nv>(C = -V2)- 



(83) 



The normalization condition (|7ip then tells us what the 
correct normalization of the GN model spinor (which dif- 
fers from the one used in the standard TDHF approach) 
should be, namely 



+ 2 - 1 - 

1 = t ct 3 </> = = "T^W- 



(84) 



Comparison with Eq. (|4"5j) shows that I = — 1/£. We now 
turn to the second spinor \ an( i the gauge transformation 
from 6*1,2 to B\ 2- The same gauge transformation fi can 
be used once again, but Eq. (|82j) is replaced by 



(85) 



A = -4C 2 



Since A on the string side is real, we have to analytically 
continue the TDHF spinors to purely imaginary spectral 
parameter ( = ±iX/2. Which sign is the correct one? 
It turns out that the normalization condition (|71l) for \ 
can only be satisfied by taking a linear combination of 
spinors with both signs [20fl . The norm vanishes if we 
keep only one sign. We choose 



X 



-i\/2)+itp(C = iX/2)] 



(86) 



where "0(C) is normalized according to (|84f before the an- 
alytic continuation. The detailed justification is given in 
the appendix. These algebraic manipulations are rooted 
in the symmetries of the GN model on the one hand and 
AdS3 space on the other hand. Neveu and Papanico- 
laou have identified a dynamical SO(2,l) symmetry of the 
classical N = 1 GN model which is also relevant for type 
I solutions of the large N quantum theory. To match the 
SO(2,l)xSO(2,l) symmetry of AdS3, one needs two in- 
dependent spinor solutions - this is the non-trivial part 
of the mapping. One of them is given directly by the 



(appropriately normalized) TDHF solution, the other in- 
volves an analytic continuation of the first one to imag- 
inary spectral parameters. This is apparently what it 
takes to embed the fermionic quantum field theory in a 
higher dimensional, classical string theory. We can now 
write down the string coordinates in compact form. In- 
troduce the basic GN TDHF spinor and the analytically 
continued one as [see Eqs. (|83I86[) ] 



0a 
06 



n~V = V(C = -a/2) 



(87) 



V2 



[VKC = -*A/2) + *V(C = *A/2)] 



normalized according to (cf. Eq. ([54")) and the appendix) 

AS 
~2~' 



Then, using the above gauge transformation, the string 
coordinates (|70|) can be expressed in the concise form 



Zi 



Ipai'a ' 



Zo 



1paH5 -0b 
i'ai'a 



(89) 



Here it looks as if the normalization of the TDHF spinors 
would simply drop out, provided we use ipa^a = 4>b4>b- 
However Eq. ([55)) for ijj a is important because it teaches 
us how to normalize the GN spinor before doing the 
analytic continuation. Otherwise, the whole procedure 
would be ill defined, since the normalization factor will 
in general depend on £. 

All of these manipulations were done under the tacit 
assumption S > 0. If S < 0, we can repeat the same 
procedure with the following changes. The gauge trans- 
formation defined in Eq. (|5Tj) has to be replaced by 



O = 



A- 1 
1X- 1 



1 



(-S)- 1 / 2 



(S < 0). 



(90) 



One can easily check that all other equations, in partic- 
ular the final result (J89J) , then remain valid. 

Using Eqs. (|88I89[) it is straightforward to verify that 
the string equations (|571 - [5^|) are satisfied. One needs the 
fact that ijja, "06 satisfy the Lax system (|52I53[) with the 
parameters h\l = A 2 /2 for ip a and h\l = — A 2 /2 for ip^; 
h^i then follows from (f55|) . The AdS3 condition (|67[) can 
be shown via the Fierz identity 

- \i> a tpb\ 2 + IVVYS'/'hT = -Tpa1pai>blpb. (91) 

The Virasoro constraints ([59^1 follow trivially from the 
useful identities 



Zi. 
Zi. 



2i 



■^li^bi, z 2 . 



X 3 S 2 

2lX I* , 7 



2A 



x 3 s 



2 ^al0bl 



62 



(92) 



Finally, the equation of motion (|68j) can be shown with 
the help of 



dY-dY 



1 

2S*2' 



(93) 
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terclockwise around its center. 



Let us illustrate this mapping between type I solutions 
of the GN model and classical strings in AdS3 by means 
of examples involving zero, one and two solitons, respec- 
tively. 

Zero soliton solution. In the vacuum, S = m = 1 
and the correctly normalized HF spinor in the present 
representation reads 




10.5 



FIG. 5: Motion of the string for the one-soliton TDHF solu- 
tion (kink at rest, A = — 2) in AdS3. The parameter t is the 
global time according to the parametrization (|96|l . The spike 
corresponds to the zero of 5* and must lie on the boundary of 
the disk. 



c 

-1/2 



i(z(-z/4() 



C=(fc-w)/2. (94) 



The spectral parameter ( is closely related to the light- 
cone momentum or energy, depending on conventions. 
The string coordinates (|89| are given by 



Zi 


= \Tie %A - cosh A+., 




z 2 


= \J~ie %A - sinhA + , 




A± 


zX z 
~2~ 2A' 


(95) 



To exhibit the motion of the string, we first introduce 
global coordinates in AdS3 through 



Z y = e" cosh £, Z 2 = sinh £ 



(96) 



where t is the global time and (£, </>) are coordinates on 
the time slice, the hyperbolic plane H2. The string at 
fixed t is most easily visualized on the Poincare disk of 
radius 1 by the choice of coordinates 



X = pcos(f>, Y = psin<j), p = J °°^ + \ - (97) 

The circumference of the disk is the boundary of AdS3 
space. Since we have two independent variables z, z while 
keeping only t fixed, Eqs. (|97|) define a one-dimensional 
curve. In the case of the vacuum, the string is evidently 
just a straight line along a diameter of the disk uniformly 
spinning around its center, see Fig. 2) This agrees with 
the findings of [l8| where the sinh-Gordon vacuum was 
mapped onto string solutions. 



One soliton solution. A kink moving with velocity v 
has the scalar potential 



S(x, t) — tanh7(.x — vt), 7 



1 



(98) 



The upper and lower components of the continuum 
spinors belonging to this potential are 



^ = M ( 1 + l Xl\ e " 



1p2 = N 



2C 
iY S 
'4? " 2C 



with 



A = kx — Lot = Qz 
2iC 



A A 



Z 



(99) 



M = - 



Y 



Y-2i(' 
1 - v 



1 + v 

Here, the normalization 

W> = -(S 



(100) 



(101) 



has been chosen in accordance with Eq. (|84|). Z\ and 
Zi can be constructed using Eqs. (|87H89|) above, but the 
result is not very transparent and will not be given here. 
In order to map out the motion of the string, we proceed 
as follows: After choosing an initial global time t, we 
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construct a trajectory in the (z, z) plane by numerically 
solving the transcendental equation 

i = arctan — — . (102) 
Y-i 

Along this trajectory we then evaluate Z\ , Z2 and convert 
them to the coordinates (X, Y) according to Eqs. (|96I97|) , 
We then plot the string and repeat the procedure for a 
sequence of time steps. We recall that we had to choose 
different gauge transformations depending on the sign of 
S, cf. Eqs. (JHIJ, (EOj). When gluing together the cor- 
responding solutions, we impose continuity on the string 
coordinates in order to specify the relative phase between 
the spinors ip a ,ipb, which would otherwise be undeter- 
mined. In this way we find that we have to change the 
sign of Z\, Z2 whenever S crosses zero. Fig. [5] shows the 
result of such a calculation, using the parameters v = 
(kink at rest) and A = — 2. As expected, the endpoints of 
the string lie on the boundary. Since S = is a singular 
point, the fold in the string also touches the boundary, 
unlike in Ref. [18j . This behavior can easily be traced 
back to the factor 1 / y/\S\ in the gauge transformation 
fl, see Eqs. ([STJ and IpTO]). Comparing Figs. [4] and [5] we 
observe that the kink string resembles the vacuum string 
at the beginning and the end of the motion shown, reflect- 
ing the characteristic behavior of the kink potential which 
also interpolates between two vacua. In order to better 
understand the shape of the strings, we also recall that 
geodesies on the Poincare disk are either straight lines 
through the center or circles intersecting the boundary 
at right angles. 

Two soliton solution. It is straightforward to repeat 
such a calculation for the kink-antikink scattering solu- 
tion discussed above, at least using Maple. Since S now 
changes sign twice, the string is obtained by patching 
together three different solutions. The resulting string 
now connects two points on the boundary of the disk 
and touches this boundary in cusps at two intermedi- 
ate points corresponding to the zeros of S. Figs. 16171 
show an example derived from kink-antikink scattering 
at v = tanh 1 w 0.76 (rapidity 1) and A = —2. Since 
the motion is already fairly complicated, we have chosen 
finer time steps than for the one soliton case. Each of 
the three parts of the string is infinitely long due to the 
Poincare metric and qualitatively follows the geodesies. 




FIG. 6: Time evolution of string derived from kink-antikink 
scattering at v = 0.76, A = —2. The representation is the 
same as in Fig. [5] but now there are two cusps on the bound- 
ary, reflecting the two zeros of S. 




FIG. 7: Continuation of Fig. [6]to larger times. 



V. SUMMARY AND DISCUSSION 

In this paper, we have reconsidered the simplest ver- 
sion of the GN model in the large N limit — the masslcss 
model with discrete chiral symmetry, Eq. (Q]). Following 
a suggestion of Dashen, Hasslacher and Neveu we have 
first studied kink-antikink scattering, using the TDHF 
language. A full analytical solution has uncovered an in- 
teresting interplay between the scalar potentials of kink 
and antikink which repell each other and the fcrmion 



constituents which travel in one direction only, being ex- 
changed during the collision. From a technical point of 
view, we have found that unlike our starting point, the 
DHN breather, the scattering solution belongs to the sim- 
ple class of type I solutions in which each single parti- 
cle orbit contributes a term proportional to the full self- 
consistent potential S(x, t) to the chiral condensate. This 
complements earlier known solutions of type I, the vac- 
uum, the kink and the kink crystal. Type I solutions 
are singled out by the fact that the TDHF equations re- 
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duce to a non-linear Dirac equation. In this way we were 
able to make contact with earlier studies of the classical 
N = 1 GN model and take advantage of the techniques 
which have been developed there. We have reproduced 
the finding of Neveu and Papanicolaou that the square of 
the self-consistent potential is related to the sinh-Gordon 
equation. The linearized version of the latter equation is 
nothing but the Klein-Gordon equation for the a meson 
of the GN model. This supports a Skyrme-type picture 
of the baryon in the non-chiral GN model. All presently 
known type I solutions can be identified (up to a triv- 
ial sign ambiguity) with well-known solitons of the sinh- 
Gordon equation, an observation which would have saved 
a lot of guesswork, had it been known before. Even more 
interesting is perhaps the fact that TV soliton solutions 
and the Lax representation of the associated linear prob- 
lem are known explicitly for the sinh-Gordon equation. 
This should enable us to generalize the two-soliton so- 
lution of the present paper to the scattering problem of 
N composite, relativistic bound states. Another possi- 
ble generalization of the present work would be to con- 
sider type II solutions, thereby covering all analytically 
known HF and TDHF solutions of both the massless and 
the massive (non-chiral) GN models. The methods of 
Ref. [17| are also developed for the N = 2 classical GN 



model, relevant for type II solutions of the large TV quan- 
tum theory. Here however, the formalism is significantly 
more involved, and one apparently needs coupled, non- 
linear differential equations for several functions, includ- 
ing the scalar potential. 

So far we have not yet mentioned the chiral GN model 
or Nambu-Jona-Lasinio model in 2 dimensions (NJL2). 
In the massless NJL2 model, transparent potentials ap- 
pear in the context of the massless baryon [24[ and the 
Shei bound state (29[. Generalizations to periodic, fi- 
nite gap potentials include the chiral spiral [301 ] and the 
twisted kink crystal |3lj . Basar and Dunne have identi- 
fied the non-linear Schrodinger equation as the relevant 
equation for the HF potential, using completely different 
techniques from the ones employed in the present paper. 
In the massive NJL2 model, approximate results from 
the derivative expansion indicate that the sine- Gordon 
equation is relevant near the chiral limit, but needs to be 
replaced by increasingly complicated, higher order differ- 
ential equations with increasing bare fermion mass 14 [. 
Here, the baryon potentials arc apparently not transpar- 
ent, so that we are dealing with neither type I nor type II 
solutions, nor with solutions of any finite type, for that 
matter. The common theme of all of these efforts is the 
search for a closed, classical theory of the self-consistent 
potential, the relevant saddle point in the functional in- 
tegral approach. In a sense, one might identify this effort 
with the search for Witten's "master field" in the large 
N limit of quantum chromodynamics [32| where unfor- 
tunetely little progress has been made until now. 

The close relationship between the GN model and the 
sinh-Gordon equation gave us the clue for yet another 
kind of mapping. It is known that the sinh-Gordon model 



and strings in AdS3 are related by some kind of gauge 
transformation, using the Pohlmcycr reduction. In the 
same vein, we have constructed a gauge transformation 
mapping type I GN solutions onto classical strings in 
AdS3. Although this has nothing to do with the cele- 
brated AdS/CFT correspondence, it is rather intriguing. 
When mapping the sinh-Gordon model to string theory, 
one has to use artificial spinor fields from the associated 
Lax pair. In our case, the spinors have a more direct 
physical meaning as TDHF solutions. Another difference 
to the string/ sinh-Gordon correspondence is a singularity 
in the gauge transformation at the zeros of S. This yields 
folded strings whose spikes must always lie on the bound- 
ary of AdS3. The existence of such a mapping seems to 
be closely related to the symmetries of the models. How- 
ever, the fact that two models have the same symmetry is 
in general not sufficient to conclude that their dynamics 
is the same. Here we have shown how to map a fermionic 
quantum field theory in the large N limit onto a classical 
string theory in a fully explicit manner. Whether this is 
of any practical use remains to be seen. If nothing else, it 
provides us with a novel way of visualizing TDHF single 
particle spinors, which contain more detailed information 
than the self-consistent potential. 
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Appendix: Definition and normalization of ipb 

Here we motivate the particular choice for the analyti- 
cally continued TDHF spinor i/jf,, Eq. (|8"T)l . Consider first 
a real spectral parameter £ . The N soliton solution of the 
TDHF equation can be factored as 



ip = tpe 



i(z(-z/4() 



(103) 



where (p is a solution of the reduced Dirac equation 

1 



2iipi. z - — (pi = Stp 2 , 
2^2,2 - 2(>2 = Sipi. 



(104) 



This merely expresses the fact that the potentials are re- 
flectionless. If ip solves Eq. (|104|) with spectral parameter 
£, then 



l5f = 



<P2 



(105) 



solves the corresponding equation with spectral parame- 
ter — £. Since we wish to introduce imaginary £ parame- 
ters, we avoid complex conjugation and write the scalar 
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condensate and the reflection property under £ — > — £ in 
the analytic form 

<Pa{-C,i) = MC-i)- (106) 



with 



We also need the normalization condition ([88)1 for xj} a 
which translates into 

$KC)v(0 = -CS. (107) 

From Eqs. (|106[) one deduces the relations 



?(-CMC) - <p(CM-0 = o, 



(108) 



In this form, all reflection properties also hold for com- 
plex £. ipi, can be written as 

i> b = -L [<p(( - ^X/2)e B + i<p(C = l X/2)e~ B ] (109) 



B 



zX 



2A' 



(110) 



Using Eqs. t|10T[) and (| 108[) for imaginary £, the conden- 
sate of f/'b becomes 



Mb = \ [^{iX/2)^(-iX/2)e 2B + <p(-i\/2)(p(i\/2)t 

+icp(iX/2)cp(iX/2) - icp(-iX/2)ip(-iX/2)] 
= icp(iX/2)ip(iX/2) 

~ 2 b > 



-IB 



(111) 



proving the 2nd half of Eq. 
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